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Abstract. We present a new efficient localized algorithm to construct, 
for any given quasi-unit disk graph G = (V, E) and any e > 0, a (1 + e)- 
spanner for G of maximum degree O(l) and total weight 0(ui(M ST)), 
where ui(MST) denotes the weight of a minimum spanning tree for V . We 
further show that similar localized techniques can be used to construct, 
for a given unit disk graph G = (V, E), a planar Cde!(l+£)(l+§)-spanner 
for G of maximum degree 0(1) and total weight 0(uj{M ST)). Here Cdei 
denotes the stretch factor of the unit Delaunay triangulation for V . Both 
constructions can be completed in O(l) communication rounds, and re- 
quire each node to know its own coordinates. 



1 Introduction 

For any fixed a, < a < 1, a graph G = (V,E) is an a- quasi unit disk graph 
(cv-QUDG) if there is an embedding of V in the Euclidean plane such that, for 
every vertex pair u, v S V, uv G E if \uv\ < a, and \uv\ ^ E if \uv\ > 1. The 
existence of edges with length in the range (a, 1] is specified by an adversary. If 
a = 1, G is called a unit disk graph (UDG). a-QUDGs have been proposed as 
models for ad-hoc wireless networks composed of homogeneous wireless nodes 
that communicate over a wireless medium without the aid of a fixed infrastruc- 
ture. Experimental studies show that the transmission range of a wireless node 
is not perfectly circular and exhibits a transitional region with highly unreliable 
links [34] (see for example Fig. la, in which the shaded region represents the 
actual transmission range). In addition, environmental conditions and physical 
obstructions adversely affect signal propagation and ultimately the transmission 
range of a wireless node. The parameter a in the a-QUDG model attempts to 
take into account such imperfections. 

Wireless nodes are often powered by batteries and have limited memory 
resources. These characteristics make it critical to compute and maintain, at 
each node, only a subset of neighbors that the node communicates with. This 
problem, referred to as topology control, seeks to adjust the transmission power 
at each node so as to maintain connectivity, reduce collisions and interference, 
and extend the battery lifetime and consequently the network lifetime. 
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Different topologies optimize different performance metrics. In this paper we 
focus on properties such as planarity, low weight, low degree, and the spanner 
property. Another important property is low interference [5, 15, 30], which we do 
not address in this paper. A graph is planar if no two edges cross each other (i.e, 
no two edges share a point other than an endpoint). Planarity is important to 
various memoryless routing algorithms [16, 4]. A graph is called low weight if its 
total edge length, defined as the sum of the lengths of all its edges, is within a 
constant factor of the total edge length of the Minimum Spanning Tree (MST). 
It was shown that the total energy consumed by sender nodes broadcasting along 
the edges of a MST is within a constant factor of the optimum [31]. Low degree 
(bounded above by a constant) at each node is also important for balancing 
out the communication overhead among the wireless nodes. If too many edges 
are eliminated from the original graph however, paths between pairs of nodes 
may become unacceptably long and offset the gain of a low degree. This renders 
necessary a stronger requirement, demanding that the reduced topology be a 
spanner. Intuitively, a structure is a spanner if it maintains short paths between 
pairs of nodes in support of fast message delivery and efficient routing. We define 
this formally below. 

Let G = (V, E) be a connected graph representing a wireless network. For 
any pair of nodes u, v € V, let sp G (u,w) denote a shortest path in G from u to 
v, and let |sp G (w, v)\ denote the length of this path. Let H C G be a connected 
subgraph of G. For fixed t > 1, H is called a t- spanner for G if, for all pairs of 
vertices u, v G V, \sp H (u,v)\ < t ■ |sp G (w,v)|. The value t is called the stretch 
factor of H . If t is constant, then H is called a length spanner, or simply a 
spanner. A triangulation of V is a Delaunay triangulation, denoted by Del(V), 
if the circumcirclc of each of its triangles is empty of nodes in V. 

Due to the limited resources and high mobility of the wireless nodes, it is 
important to efficiently construct and maintain a spanner in a localized manner. 
A localized algorithm is a distributed algorithm in which each node u selects all 
its incident edges based on the information from nodes within a constant number 
of hops from u. Our communication model is the standard synchronous message 
passing model, which ignores channel access and collision issues. In this commu- 
nication model, time is divided into rounds. In a round, a node is able to receive 
all messages sent in the previous round, execute local computations, and send 
messages to neighbors. We measure the communication cost of our algorithms in 
terms of rounds of communication. The length of messages exchanged between 
nodes is logarithmic in the number of nodes. 

Our Results. In this paper we present the first localized method to construct, 
for any QUDG G = (V, E) and any e > 0, a (1 + e)-spanncr for G of maximum 
degree 0(1) and total weight 0(uj{MST)), where uj{MST) denotes the weight 
of a minimum spanning tree for V. We further extend our method to construct, 
for any UDG G — (V,E), a planar spanner for G of maximum degree 0(1) and 
total weight 0(lu(MST)). The stretch factor of the spanner is bounded above by 
Cdez(L+ £ )(l. + § ), where Cdei is the stretch factor of the unit Delaunay triangula- 
tion for V (Cdei < 2.42 [20]). This second result resolves an open question posed 



by Li et al. in [22]. Both constructions can be completed in 0(1) communication 
rounds, and require each node to know its own coordinates. 



1.1 Related Work 

Several excellent surveys on spanners exist [27,26,14,25]. In this section we 
restrict our attention to localized methods for constructing spanners for a given 
graph G = (V,E). We proceed with a discussion on non-planar structures for 
UDGs first. Existing results are summarized in the first four rows of Table 1. 

The Yao graph [33] with an integer parameter k > 6, denoted YGk, is defined 
as follows. At each node u £ V, any k equal-separated rays originated at u define 
k cones. In each cone, pick a shortest edge uv, if there is any, and add to YGk the 
directed edge uv. Ties are broken arbitrarily or by smallest ID. The Yao graph is a 
spanner with stretch factor 1 _ 2s h 17r /fc > however its degree can be as high as n— 1. 
To overcome this shortcoming, Li et al. [18] proposed another structure called 
YaoYao graph YYk, which is constructed by applying a reverse Yao structure 
on YGk- at each node u in YGk, discard all directed edges vu from each cone 
centered at u, except for a shortest one (again, ties can be broken arbitrarily 
or by smallest ID). YYk has maximum node degree 2k, a constant. However, 
the tradeoff is unclear in that the question of whether YYk is a spanner or not 
remains open. Both YGk and YYk have total weight 0(n) ■ ui(MST) [6]. Li 
et al. [32] further proposed another sparse structure, called YaoSink YSk, that 
satisfies both the spanner and the bounded degree properties. The sink technique 
replaces each directed star in the Yao graph consisting of all links directed into a 
node u, by a tree T[u) with sink u of bounded degree. However, neither of these 
structures has low weight. 
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Table 1. Results on localized methods for UDGs. 



We now turn to discuss planar structures for UDGs. The relative neighbor- 
hood graph (RNG) [29] and the Gabriel graph (GG) [12] can both be constructed 
locally, however neither is a spanner [2]. On the other hand, the Delaunay tri- 
angulation Del(V) is a planar t-spanner of the complete Euclidean graph with 
vertex set V. This result was first proved by Dobkin, Friedman and Supowit [11], 
for t = 1+ 2 v ^ 7r w 5.08, and was further improved to t = ^p7r ~ 2.42 by Keil 
and Gutwin [17]. Das and Joseph [7] generalize these results by identifying two 
properties of planar graphs, the good polygon and diamond properties, which 
imply that the stretch factor is bounded above by a constant. 

For a given point set V, the unit Delaunay triangulation of V , denoted 
UDel(V), is the graph obtained by removing all Delaunay edges from Del(V) 
that are longer than one unit. It was shown that UDel(V) is a t-spanner of the 
unit-disk graph UDG0O, with t = ^tt w 2.42 [20]. 

Gao et al. [13] present a localized algorithm to build a planar spanner called 
restricted Delaunay graph (RDG), which is a supergraph of UDel(V). Li et al. [20] 
introduce the notion of a k-localized Delaunay triangle: Aabc is called k-localized 
Delaunay if the interior of its circumcircle does not contain any node in V 
that is a fc-neighbor of a, b or c, and all edges of Aabc are no longer than one 
unit. The authors describe a localized method to construct, for fixed k > 1, 
the k- localized Delaunay graph LDel k (V), which contains all Gabriel edges and 
edges of all k- localized Delaunay triangles. They show that (i) LDel k (V) is a 
supergraph of UDel(V) (and therefore a ^p7r-spanncr), (ii) LDel fe (V) is pla- 
nar, for any k > 2, and (iii) LDel 1 (F) may not be planar, but a planar subgraph 
PLDel(F) C LDel 1 (V) that retains the spanner property can be locally extracted 
from LDel 1 (T^). Their planar spanner constructions take 4 rounds of communi- 
cation and a total of 0(n) messages (O(nlogn) bits). Araujo and Rodrigues [1] 
improve upon the communication time for PLDel and devise a method to com- 
pute PLDel(V) in one single communication step. Both PLDel(T^) and LDel fc (y), 
for k > 1, may have arbitrarily large degree and weight. 

To bound the degree, several methods apply the ordered Yao structure on 
top of an unbounded-degree planar structure. This idea was first introduced by 
Bosc et al. in [3], and later refined by Li and Wang in [32, 23]. Since the ordered 
Yao structure is relevant to our work in this paper as well, we pause to discuss 
the Ordered Yao method for constructing this structure. The Ordered Yao 
method is outlined in Table 2. The main idea is to define an ordering ir of the 
nodes such that each node u has a limited number of neighbors (at most 5) who 
are predecessors in 7r; these predecessors are used to define a small number of 
open cones centered at u, each of which will contain at most one neighbor of u 
in the final structure. To maintain the spanner property of the original graph, 
a short path connecting all neighbors of u in each cone is used to replace the 
edges incident to u that get discarded from the original graph. 

Thm. 1 summarizes the important properties of the structure computed by the 
Ordered Yao method. 



Theorem 1. If G is a planar graph, then the output G' obtained by executing 
Ordered Yao(G) is a planar (1 + ■§)- spanner for G of maximum degree 25 [32]. 



Algorithm Ordered Yao(G = (V,E)) [32] 



{1. Find an order n for V:} 

Initialize i = 1 and G; = G. 
Repeat for i = 1, 2, . . . , \V\ 

Remove from G; the node u of smallest degree 

(break ties by smallest ID.) 

Call the remaining graph Gj+i. 

Set ty u = n — i + 1. 

{2. Construct a bounded-degree structure for G:} 

Mark all nodes in V unprocessed. Initialize E' <— and G' = (V, -E'). 
Repeat \V\ times 

Let u be the unprocessed node with the smallest order ir u . 
Let vi, V2, ■ ■ ■ , Vh be the be the processed neighbors of u in G (h < 5). 
Shoot rays from u through each Vi, to define /i sectors centered at u. 
Divide each sector into fewest open cones of degree at most it/3. 
For each such open cone C u (refer to Fig. above) 

Let si, S2, ■ ■ ■ , s m be the geometrically ordered neighbors of u in C u . 

Add to E' the shortest us; edge. 

Add to E' all edges SjSj+i, for j = 1, 2, . . . , m — 1. 
Mark node u processed. 

Output G = (V,E r ). 



Table 2. The Ordered Yao method. 




Song et al. [28] apply the ordered Yao structure on top of the Gabriel graph 
GG(V) to produce a planar bounded-degree structure OrdYaoGG. Their result 
improves upon the earlier localized structure YaoGG [18], which may not have 
bounded degree. Both YaoGG and OrdYaoGG are power spanners, however neither 
is a length spanner. 

The first efficient localized method to construct a bounded-degree planar 
spanner was proposed by Li and Wang in [32,23]. Their method applies the 
ordered Yao structure on top of LDel(V) to bound the node degree. The resulted 
structure, called BPS(V) (Bounded-Degree Planar Spanner), has degree bounded 
above by 19 + [^] , where < a < | is an adjustable parameter. The total 
communication complexity for constructing BPS(F) is 0(n) messages, however 
it may take as many as 0{n) rounds of communication for a node to find its rank 



in the ordering of V (a trivial example would be n nodes lined up in increasing 
order by their ID). The BPS structure does not have low weight [19]. 

The first localized low-weight planar structure was proposed in [19]. This 
structure, called RNG' , is based on a modified relative neighborhood graph, and 
satisfies the planarity, bounded-degree and bounded-weight properties. A similar 
result has been obtained by Li, Wang and Song [22], who propose a family of 
structures, called Localized Minimum Spanning Trees LMSTfc, for k > 1. The 
authors show that LMSTfc is planar, has maximum degree 6 and total weight 
within a constant factor of u>(MST), for k > 2. Their result extends an earlier 
result by Li, Hou and Sha [24], who propose a localized MST-based method to 
compute a local minimum spanning tree structure. However, neither of these 
low- weight structures satisfies the spanner property. Constructing low- weight, 
low-degree planar spanners in few rounds of communication is one of the open 
problems we resolve in this paper. 



2 Our Work 



We start with a few definitions and notation to be used through the rest of the 
paper. For any nodes u and v, let uv denote the edge with endpoints u and v; 
uv is the edge directed from u to v; and \uv\ denotes the Euclidean distance 
between u and v. Let C u denote an arbitrary cone with apex u, and let C u (v) 
denote the cone with apex u containing v. For any edge set E and any cone C u , 
let E fl C u denote the subset of edges in E incident to u that lie in C u . 

We assume that each node u has a unique identifier ID(tt) and knows its 
coordinates (x u ,y u )- Define the identifier ID(mu) of a directed edge uv to be 

the triplet (\uv\, ID(tt), ID(u)). For any pair of directed edges uv and u'v', we 
say that ID(uv) < ID(wV) if and only if one of the following conditions holds: 
(1) \uv\ < \u'v'\, or (2) \uv\ = \u'v'\ and ID(u) < ID(u'), or (3) \uv\ = \u'v'\ 
and ID(m) = ID(u') and ID(u) < ID(u'). For an undirected edge uv, define 
IT)(uv) — min{lD(wu), ID(vu)}. Note that according to this definition, each edge 
has a unique identifier. 

Let H — (V, Eh) be an arbitrary subgraph of G = (V, E). A subset L u <zV 
is an r-cluster in H with center u if, for any v e L u , \sp H (u, v)\ < r. A set 
of disjoint r-clusters {L Ul ,L U2 , . . .} form an r-cluster cover for V in H if they 
satisfy two properties: (i) for i ^ j, \sp H (ui,Uj)\ > r (the r-packing property), 
and (ii) the union UiL Ui covers V (the r-covering property). 

For any node subset U C V, let G[U] denote the subgraph of G induced by 
U. A set of node subsets V\, V2, . ■ ■ C V is a clique cover for V if the subgraph 
of G[Vi] is a clique for each i, and U^ =1 Fi = V. 

The aspect ratio of an edge set E is the ratio of the length of a longest edge 
in E to the length of a shortest edge in E. The aspect ratio of a graph is defined 
as the aspect ratio of its edge set. 



2.1 The LOS Algorithm 



In this section we describe an algorithm called LOS (Localized Optimal Spanner) 
that takes as input an a-QUDG G = (V, E), for fixed < a < 1, and a value 
e > 0, and computes a (1 +e)-spanner for G of maximum degree O(l) and total 
weight 0(uj(MST)). The main idea of our algorithm is to compute a particular 
clique cover Vi, V2, ■ ■ ■ for V, construct a (1 + e)-spanner for each G[Vi], then 
connect these smaller spanners into a (1 + e)-spanner for G using selected Yao 
edges. In the following we discuss the details of our algorithm. 
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Fig. 1. (a) The a-QUDG model (b) Constructing a clique cover for V (c) Clique 
ordering. 



Let < < -j= and < S < /3/4 be small constants to be fixed later. 
To compute a clique cover for V , we start by covering the plane with a grid of 
overlapping square cells of size (3 x /3, such that the distance between centers of 
adjacent cells is — 28. Note that any two adjacent cells define a small band of 
width S where they overlap. The reason for enforcing this overlap is to ensure 
that edges not entirely contained within a single grid cell are longer than 5, i.e., 
they cannot be arbitrarily small. We identify each grid cell by the coordinates 
(i,j) of its upper left corner. Any two vertices that lie within the same grid cell 
are no more than a distance apart and therefore are connected by an edge in 
G. This implies that the collection of vertices in each non-empty grid cell can 
be used to define a clique element of the clique cover. We call this particular 
clique cover a ((3 , S)- clique cover. Let Vi,V2, ■ ■ ■ be the elements of the ((3,5)- 
clique cover for V. Note that, since 5 < f3/4, a node u can belong to at most 
four subsets Vi. 

Our LOS method consists of 4 steps. First we construct, for each G[Vi], a 
(1 + e)-spanner of degree O(l) and weight 0(u(MST(Vi)). Various methods for 
constructing Hi exist - for instance, the well-known sequential greedy method 
produces a spanner with the desired properties [8]. Second, we use the Yao 
method to generate (1 + e) -spanner paths between longer edges that span differ- 
ent grid cells. Third, we apply the reverse Yao step to reduce the number of Yao 
edges incident to each node. Finally, we apply a filtering method to eliminate 
all but a constant number of edges incident to a grid cell. This fourth step is 
necessary to ensure that the output spanner has bounded weight. These steps 



Algorithm LOS(G = (V,E),e) 

{1. Compute a (1 + e)-spanner cover:} 

Fix < (3 < -z= and < <5 < (3/4. 

Compute a (/3, (5)-clique cover Vi, V2, . . . for V. 

For each i, compute a (1 + e)-spanner Hi for G[Vi] using the method from [8]. 
Initialize H = Uji/j. Let Eo = {uv £ E \ uv £ G[Vi] for any i}. 

{2. Apply Yao on E< } :} 

Let k be the smallest integer satisfying cos 2zl — sin 2zl > _ 
For each node w, divide the plane into k incident equal-size cones. 
Initialize Ey <— 0. 

For each cone C u such that i?o fl C u is non-empty 

Pick the edge uv G E PI C u of smallest ID and add uv to £V . 

{3. Apply reverse Yao on Ey:} 
Initialize Eyy <— _Ey. 

For each cone C u such that Ey n C u is non-empty 

Discard from _Ey all edges vu G -Ey fl C u , but the one of smallest ID. 

{4. Select connecting edges from Eyy:} 

Pick r such that r < g±lKl±£ii£g« fl-Binfl)-(g+i+ e ) > where fl = 2?r/ ^ 

Compute an r-cluster cover for V in ff . 

Let Ei C £yy contain all Yao edges connecting cluster centers. Add E\ to _ff. 
Output H = (V,E H ). 



Table 3. The LOS algorithm. 



are described in detail in Table 3. Note that the Yao and reverse Yao steps are 
restricted to edges in the set E whose aspect ratio is bounded above by 1/5. 
The next three theorems prove the main properties of the LOS algorithm. 

Theorem 2. The output H generated by LOSfG, e) is a (1 + e)- spanner for G. 

Proof. Let uv € E be arbitrary. If uv e G[Vi] for some i, then Hi C H contains 
a (1 + e)-spanner uw-path (since Hi is a (1 + e)-spanner for G[Vi]). Otherwise, 
uv e E . The proof that H contains a (1 + e)-spanner uw-path is by induction on 
the ID of edges in E . Let uv e E n be the edge with the smallest ID and assume 
without loss of generality that TD(uv) = TD{uv). Since ID(uv) is smallest, uv 
gets added to Ey in step 2, and it stays in Eyy in step 3. If uv <G H at the end 
of step 4, then sp H (u, v) — uv. Otherwise, let ab be the edge selected in step 4 
of the algorithm, such that u G L a and v € L& (see Fig. 2a). Since L a and Lb are 
both r-clusters, we have that \sp H (u, a)\ < r and \sj> H (v,b)\ < r. It follows that 
\ua\ < r and \vb\ < r. By the triangle inequality \ab\ < \uv\ + 2r and therefore 



sp ff (u, v) < \ab\ + 2r < \uv\ + Ar < (1 + e)|w|, for any r < Se/4 (satisfied by 
the r values restricted by the algorithm). This concludes the base case. 

To prove the inductive step, let uv € Eq be arbitrary, and assume that H 
contains (1 + e)-spanner paths between the endpoints of any edge whose ID is 
lower than ID(uv). 




(a) (b) 

Fig. 2. Thm. 2: (a) Base case, (b) sp H (it, ui) ©sp H (iti, a) ©afc©sp /J -(6, Vi) © sj> H (vi, v) 
is a (1 + e)-spanner wv-path. 



Let uvi G C u (v) be the Yao edge selected in step 2 of the algorithm; let 
uivi € C Vl (u) be the Yao Yao edge selected in step 3 of the algorithm; and let 
ab 6 H be the edge added to H in step 4 of the algorithm, such that ui € L a 
and vi € Lb (see Fig. 2b). Note that u and m may be disjoint or may coincide, 
and similarly for v and v±. In either case, the chains of inequalities ID(ui^i) < 
ID(wi) < IV (uv) and < < \uv\ hold. Let u[ be the projection of u\ 

on uv\. By the triangle inequality, 

\uu\\ < \uu[\ + \u[ui\ — \uvi\ — \uiVi\ + \uiUi\ < \uvi\ — \uivi\co$6+\uivi\siri6. 

(1) 

Similarly, if v[ is the projection of v\ on uv, we have 

\v\v\ < \vv[\ + \v[vx\ = \uv\ — \uv[\ + \v[vi\ < \uv\ — \uvi\ cos + \uvi\ sin9. (2) 

Since \uu\\ < \uvi\ < \uv\ and \viv\ < \uv\, by the inductive hypothesis H 
contains (l + e)-spanner paths sp H (u, Ui) and sp H (vi,v). Let Pj = sp H (u, Ui)(B 
sp H (vi, v). The length of Pi is 

\Pi\ <(l + e)-(\uu 1 \ + \v 1 v\). 

Substituting inequalities (1) and (2) yields 

\Pi\ < (l + £)\uv\ + (l + £)\uv 1 \(l-cos6 + sm0)-(l+e)\uiv 1 \(cos6-sm6). (3) 

Next we show that the path P = P\® sj> H (ui,a) (Bab® sp H (fe, V\) is a (1 + e)- 
spanner path from u to v in H, thus proving the inductive step. Using the fact 
that \ab\ < 2r + \uiV\\, \sj> H (u±, a)\ < r and \sp H (b, vi)\ < r, we get 



\P\ < |Pi| + K«x|+4r. 



(4) 



Substituting further > 5 and \uv\\ < 1 in (3) and (4) yields 

|P| < (1 + e)\uv\ + (4r + (1 + s)(l - cos 6 + sin 6) - 5(1 + s) (cos 6 - sin 6) - 5). 

Note that the second term on the right side of the inequality above is non-positive 
for any r and satisfying the conditions of the algorithm: 



Before proving the other two properties of H (bounded degree and bounded 
weight), we introduce an intermediate lemma. For fixed c > 0, call an edge 
set F c-isolated if, for each node u incident to an edge e e F, the closed disk 
disk^c) centered at u of radius c contains no other endpoints of edges in F. 
This definition is a variant of the isolation property introduced in [10]. Das et al. 
show that, if an edge set F satisfies the isolation property, then u)(F) is within 
a constant factor of the minimum spanning tree connecting the endpoints of F. 
Here we prove a similar result. 

Lemma 1. Let F be a c-isolated set of edges no longer than 1. Then u)(F) = 
0(1) ■ w(T), where T is the minimum spanning tree connecting the endpoints of 
edges in F. 

Proof. Let P be a Hamiltonian path obtained by a taking a preorder traversal 
of T. If each edge uv e P gets associated a weight value u(uv) = |sp T (w, v)\, 
then it is well-known that uj(P) < 2ui(T). So in order to prove that w(F) is 
within a constant factor of w(T), it suffices to show that uj(F) = 0(lo(P)). Since 
F is c-isolated, the distance between any two vertices in T is greater than c and 
therefore w(P) > (n — l)c. On the other hand, no edge in F is greater than 1 
and therefore uo(F) < n. It follows that uu(F) = 0(w(P)). □ 

Theorem 3. The output H generated by running LOS fG, t) has maximum de- 
gree 0(1) and total weight 0(1) • lu(MST). 

Proof. The fact that H has maximum degree O(l) follows immediately from 
three observations: (a) each spanner Hi constructed in step 1 of the algorithm 
has degree O(l) [8], (b) a node u belongs to at most four subgraphs Hi, and (c) 
a node u is incident to a constant number of Yao edges (at most 2k) [18]. 

We now prove that the total weight for H is within a constant factor of 
ui(MST), which is optimal. The main idea is to partition the edge set Eh into 
a constant number of subsets, each of which has low weight. Consider first the 
(1 + e)-spanners constructed in step 1 of the algorithm. Each (1 + £r)-spanner Hi 
corresponds to a grid cell (i,j). Let F denote the set of edges in UeHe. Define the 
edge set F s C F to contain all spanner edges corresponding to those grid cells 
(i, j) whose indices i and j satisfy the condition (i mod 3) x 3 + j mod 3 = s. 
Intuitively, if two edges e\, e 2 € F s lie in different grid cells, then those grid cells 




This completes the proof. 



□ 



are separated by at least two other grid cells (see Fig. lc). This further implies 
that the closest endpoints of e\ and e-i are distance a or more apart. Also notice 
that it takes only 9 subsets F\, F2, . . . , Fg to cover F. 

Next we show that lj(F s ) = 0(u(T s )) for each s = 1, 2, . . . , 9, where T s is a 
minimum spanning tree connecting the endpoints in F s . To see this, first observe 
that F s combines the edges of several low-weight (1 + e)-spanners H Sl , H S2 , . . . 
with the property that u>(H Sl ) = 0(u>(T Sl )), where T 8l is a minimum spanning 
tree connecting the nodes in H Sl . Thus, in order to prove that ui(F s ) = 0(u(T 8 )), 
it suffices to show J2i = 0(w(T s )). We will in fact prove that 

£>(T Si ) < w(T a ) 

i 

We prove this by showing that, if Prim's algorithm is employed in constructing T s 
and T Si , then T Si C T s , for each i. Since the trees T Si are all disjoint (separated by 
at least 2 grid cells), the claim follows. Recall that Prim's algorithm processes 
edges by increasing length and adds them to T s as long as they do not close 
a cycle. This means that all edges shorter than a are processed before edges 
longer than a. Let e G T Si be arbitrary. Then |e| < a, since T Si is restricted 
to one grid cell only of diameter a. If e ^ T s , then it must be that e closes a 
cycle C at the time it gets processed. Note however that C must lie entirely in 
the grid cell containing T Si , since C contains edges no longer than a, and all 
edges with endpoints in different cells are longer than a. Furthermore, C must 
contain an edge e' ^ T Si such that |e'| < |e|. The case |e'| = \e\ cannot happen 
if Prim breaks ties in the same manner in both T s and T Si , so it must be that 
|e'| < \e\. But then we could replace e in T Si by e', resulting in a smaller spanning 
tree, a contradiction. It follows that e € T s and therefore T Si C T s , for each i. 
This concludes the proof that lj(F s ) = 0(oj(T s )) 7 for each s. Since there are 
at most 9 such sets F s that cover F and since oj(T s ) < uj(MST), we get that 
Lo(F) = 0(lo(MST)). 

It remains to prove that u{E H \ F) = 0{lo{MST)). Let d < 2k be the 
maximum number of edges in Eh \ F incident to any node in H. Partition 
the edge set Eh \ F into no more than 2d < Ak subsets Ei 7 E 2} ■ ■ ., such that 
no two edges in Ei share a vertex, for each i. We now show that co(Ei) = 
0(uj(MST)), for each i. Since there are only a constant number of sets Ei (4fc 
at most), it follows that uj{E h \ F) = 0(lj(MST)). The key observation to 
proving that iv(Ei) = 0{uj(MST)) is that any two edges uv, ab € E t have their 
closest endpoints - say, u and a - separated by a distance of at least r/t. This 
is because t\ua\ > \sp H (u, a)\ > r; the first part of this inequality follows from 
the spanner property of H, and the second part follows from the fact that u 
and a are centers of different r-clusters (a property ensured by step 4 of the 
algorithm). This implies that Ei is r/t- isolated, and by Lem. 1 we have that 
u{Ei) = 0(uj{MST)). 

We have established that u){F) = 0{uo{MST)) &nduo{E H \F) = 0(lu(MST)). 
It follows that w(H) = w(E H ) = 0(uj(MST)) and this completes the proof. □ 



Theorem 4. The LOS algorithm can be implemented in 0(1) rounds of com- 
munication using messages that are O(logn) bits each. 

Proof. Let x u and y u denote the coordinates of a node u. At the beginning of the 
algorithm, each node u broadcasts the information (lD(u), x u , y u ) to its neighbors 
and collects similar information from its neighbors. Each node u determines the 
grid cell(s) (i, j) it belongs to from two conditions, ia/^/2 < x u < (i + l)a/y/2 
and ja/\/2 < y u < (j + l)a/\/2. Similarly, for each neighbor v of u, each node 
u determines the grid cell(s) that v belongs to. Thus step 1 of the algorithm can 
be implemented in one round of communication: using the information from its 
neighbors, each node u computes the clique corresponding to those cells 
that u belongs to (at most 4 of them), then u computes a (1 + £)-spanner for 
each such clique by performing local computations. Note that knowledge of node 
coordinates is critical to implementing step 1 efficiently. 

Step 2 (the Yao step) and step 3 (the reverse Yao step) of the algorithm are 
inherently local: each node u computes its incident Yao and Yao Yao edges based 
on the information gathered from its neighbors in step 1. 

It remains to show that step 4 can also be implemented in 0(1) rounds of 
communication. We will in fact show that eight rounds of communication suffice 
to compute an r-cluster cover for V in H . Define U s to be the set of vertices 
that lie in the grid cells such that (i mod 2) x 2 + j mod 2 = s. This is 
the same as saying that two vertices that lie in different cells are about one grid 
cell apart. Note that V — Uf =1 ?7 s . To compute an r-cluster cover for V, each 
node u executes the ClusterCover method described below. For simplicity we 
assume that r > 5, so that two cluster centers that lie in different grid cells are 
at least distance r apart. However, the ClusterCover method can be easily 
extended to handle the situation r < S as well. 



Computing a ClusterCover^, r) 

Repeat for s = 1, 2, 3, 4 

(A) Collect information on cluster centers from neighbors (if any). 
If u belongs to U s 

Let Ve C Us be the clique containing u (computed in step 1 of LOS). 

(B) Broadcast information on existing cluster centers in Ve to all nodes in V(. 

(C) For each existing cluster center w € Ve 

Add to C m all uncovered nodes v £ Ve such that sp H (w,v) < r. 
Mark all nodes in C w covered. 

(D) While Ve contains uncovered nodes 

Pick the uncovered node w € Ve of highest ID. 

Add to C w all uncovered nodes v £ Ve such that sTp H (v,w) < r. 

Mark all nodes in C w covered. 

(E) Broadcast the cluster centers computed in step (C) to all neighbors. 



No information on existing cluster centers is available in the first iteration of the 
ClusterCover method (i.e, for s = 1). Each node in Ui skips directly to step 



(D), which implements the standard greedy method for computed an r-clique 
cover for a given node set (Vg in our case). In the second iteration, some of the 
clusters computed during the first iteration might be able to grow to incorporate 
new vertices from L^- This is particularly true for cluster centers that lie in the 
overlap area of two neighboring cells. Information on such cluster centers is 
distributed to all relevant nodes in step (E) in the first iteration, then collected 
in step (A) and forwarded to all nodes in Vi in step (B) in the second iteration. 
This guarantees that all nodes in have a consistent view of existing cluster 
centers in Ve at the beginning of step (C). Existing clusters grow in step (C), if 
possible, and new clusters get created in step (D), if necessary. This procedure 
shows that it takes no more than 8 rounds of communication to implement step 
4 of the LOS algorithm. One final note is that information on a constant number 
of cluster centers is communicated among neighbors in steps (A), (B) and (D) 
of the ClusterCover method. This is because only a constant number of r- 
clusters can be packed into a grid cell. So each message is O(logn) bits long, 
necessarily so to include a constant number of node identifiers, each of which 
takes O(logn) bits. □ 



2.2 The PLOS Algorithm 

In this section we impose our spanner to be planar, at the expense of a bigger 
stretch factor. This tradeoff is unavoidable, since there are UDGs that contain 
no (1 + e)-spanner planar subgraphs, for arbitrarily small e (a simple example 
would be a square of unit diameter). 

Our PLOS algorithm consists of 4 steps. In a first step we construct the unit 
Delaunay triangulation UDel(V) using the method described in [21]. Remaining 
steps use the grid-based idea from Sec. 2.1 to refine the Delaunay structure. 
Let Vi,V2,... be a (/3, <S)-clique cover for V, as defined in Sec. 2.1. In step 2 
of the algorithm we apply the Ordered Yao method on edge subsets of UDel 
incident to each clique Vi . The reason for restricting this method to each clique, 
as opposed to the entire spanner Wel(V) as in [32], is to reduce the total of 0(n) 
rounds of communication to O(l). The individual degree of each node increases 
as a result of this alteration, however it remains bounded above by a constant. 
Steps 3 and 4 aim to reduce the total weight of the spanner. Step 3 uses a Greedy 
method to filer out edges with both endpoints in one same clique Vi . Step 4 uses 
clustering to filter out edges spanning multiple cliques. These steps are described 
in detail in Table 4. The reason for breaking up step 3 of the algorithm into 4 
different rounds (for k = 1, ... ,4) will become clear later, in our discussion of 
communication complexity (Thm. 9). We now turn to proving some important 
properties of the output spanner. We start with a preliminary lemma. 

Lemma 2. The graph YDel constructed in step 2 of the PLOS algorithm is a 
planar t\- spanner for G, for any t\ > G'deii^ + 1)- Furthermore, for each edge 
ah e G, YDel contains a t\-spanner ab-path with all edges shorter than ab [32]. 



Algorithm PLDS(G = (V, E),e) 

{1. Start with the localized Delaunay structure for G:} 
Compute LDel = (V,E LBel ) for G using the method from [21]. 
Fix < (3 < -±= and < <5 < f . Compute a (/?, <5)-clique cover Vi, Vi, . . . for V. 

{2. Bound the degree:} 

For each clique Vi do the following: 

2.1 Let Ei C -BuDei contain all unit Delaunay edges incident to nodes in Vi. 

2.2 Execute YDeli <- Ordered YAo(Gi = (V,Ei)) (see Table 2). 
Set YDel = (V, E Wel ) = \J f YDel,. 

{3. Bound the weight of edges confined to single grid cells:} 

Initialize Eh = and H = (V, Eh). 
Repeat for k = 1, 2, 3,4 

{Use Greedy on non-adjacent grid cells:} 

For each grid cell L = L(i,j) such that (i mod 2) x 2 + j mod 2 = k 

3.1 Let = £rDei H L contain all edges in YDel that lie in L. 

Let Eq — E mel \ El and Q — (V, Eq) define the query graph for El- 

3.2 Sort El in increasing order by edge ID. 

For each edge e = uv £ El, resolve a shortest path query: 
If spg(u,w) > (1 +e)|wu| then add uv to H and Q. 
Otherwise, eliminate uv from YDel. 

{4. Bound the weight of edges spanning multiple grid cells:} 

Pick r such that r < ^ and compute an r-cluster cover for YDel. 
Add to H those edges in YDel connecting cluster centers. 

Output H = (V,E H ). 



Table 4. The PLOS algorithm. 



Proof. LDel is a planar C^-spanncr for G [21]. By Thm. 1, YDel^ is a planar 
(| + l)-spanner for Gi, for each i. These together with the fact that LDel = \J t d 
show that YDel is a ii-spanner for G. 




(a) c w (b) 



Fig. 3. YDel is planar: edges ab and uv cannot cross. 

The fact that YDel is planar follows an observation in [32] stating that, if 
a non-Delaunay edge e e YDel crosses a Delaunay edge e', then e' must be 
longer than one unit and does not belong to YDel. More precisely the following 
properties hold: 

(a) A non-Delaunay edge ab € YDel cannot cross a Delaunay edge uv € YDel. 
Recall that each non-Delaunay edge ab € YDel closes an empty triangle Aabc 
whose other two edges ac and be are Delaunay edges. Thus, if ab crosses uv, 
then at least one of ac and be must cross uv, contradicting the planarity of 
LDel (see Fig 3a). 

(a) No two non-Delaunay edges ab,uv G YDel cross. The arguments here are 
similar to the ones above: if ab and uv intersect, then at least two of the 
incident Delaunay edges intersect, contradicting the planarity of LDel (see 
Fig. 3b). 

The second part of the lemma follows from [32]. □ 

Theorem 5. The output H generated by PLOS (G,e) is a planar t- spanner for 
G, for any constant t > Cd e ;(l + e)(l + §). 

Proof. Since H C YDel, by Lem. 2 we have that H is planar. We now show that 
H is a i-spanner for G. The proof is by induction on the length of edges in H . 
The base case corresponds to the edge uv e G of smallest ID. Clearly uv e LDel, 
since uv is a Gabriel edge. Also uv G YDel, since it has the smallest ID among all 
edges and therefore it belongs to the Yao structure for LDel. We now distinguish 
two cases: 

(a) There is a grid cell containing both u and v. In this case uv G H , since uv is 
the first edge queried by Greedy in step 3 and therefore it gets added to H. 

(b) There is no grid cell containing both u and v. Let ab be the edge selected 
in step 4 of the algorithm, such that m£L„ and v € Lb (see Fig. 2a). Then 
arguments similar to the ones used for the base case of Thm. 2 show that 
sp H (u, a) (Bab® sp H (b,v) is a (1 + e)-spanner uv-p&th, for any r < £(5/4. 



This concludes the base case. To prove the inductive step, pick an arbitrary edge 
uv G G, and assume that H contains i-spanner paths between the endpoints of 
each edge in G of smaller ID. By Lem. 2, YDel contains a y^-spanner path 
u = w , Mi, ...,u s = v: 

s t 

< rqr^lH ( 5 ) 

i=0 

For each edge UiU i+ i e YDel, one of the following cases applies: 

(a) There is a grid cell containing both m and Ui+\. In this case, the Greedy step 
(step 3 of the algorithm) guarantees that |sp ff (m, < (1 + e)\uiUi + i\. 

(b) There is no grid cell containing both Ui and u i+ \. Arguments similar to the 
ones for the base case show that |sp ff (m, Ui + i)\ < (1 + e)\uiU i+ i\. 

In either case, H contains a (1 + e)-spanner UiU i+ i-path. This together with (5) 
shows that 

s s 

|sp H (u, v)\ = ^ \ s Ph( u ^ u *+i)\ < (1 + e) X] ^ 

i=0 i=0 

This completes the proof. □ 

Theorem 6. The output H generated by PLDS has maximum degree 0(1). 

Proof. Since H C YDel, it suffices to show that the graph YDel constructed 
in step 2 of the PLOS algorithm has degree bounded above by a constant. By 
Thm. 1, the maximum degree of YDel^ is 25, for each %. Also note that unit disk 
centered at a node u intersects 0(^2) grid cells, meaning that u is a neighbor of 
nodes in O(-^j-) grid cells and therefore belongs to a constant number of graphs 
YDel,;. This implies that the maximum degree of u is 25 • O(^j-), which is a 
constant. □ 

Definition 1. [Leapfrog Property] For any t > t' > 1, a set F of edges has 
the (t', t)-leapfrog property if, for every subset S = {uiVi,u 2 V2, ■ ■ ■ , u m v m } of F, 

t' ■ \uivi\ < X \uiVi\ +t ■ ( + \v m ui\y (6) 

j=2 i=l 

Das and Narasimhan [9] show the following connection between the leapfrog 
property and the weight of the spanner. 

Lemma 3. Let t > t' > 1. If the line segments F in d-dimensional space 
satisfy the (t' ,t) -leapfrog property, then lu(F) = 0(lu(MST)), where MST is a 
minimum spanning tree connecting the endpoints of line segments in F. 

Lemma 4. At the end of each iteration k in step 3 of the PLOS algorithm, for 
k = 1, . . . ,4, Q contains (1 + e) k -spanner paths between the endpoints of any 
YDel edge processed in iterations 1 through k. 



Proof. The proof is by induction on k. The base case corresponds to k = 1. In 
this case, Greedy ensures that Q contains a (l + e)-spanner ww-path for each edge 
uv processed in this iteration. This is because uv <G YDel either gets added to Q 
in step 3.1 (and never removed thereafter), or gets queried in step 3.2. To prove 
the inductive step, consider a particular iteration k > 1, and assume that the 
lemma holds for iterations I = 1 ... k — 1. Again Greedy ensures that Q contains 
a (l+e)-spanner iru-path for each edge uv processed in iteration k. Consider now 
an arbitrary edge uv processed in iteration £ < k. By the inductive hypothesis, 
at the end of round fc— 1, Q contains a (1 + e) fc_1 -spanner path p(u, v). However, 
it is possible that p(u, v) contains edges processed in round k (since Greedy does 
not restrict p{u, v) to lie entirely in the cell containing uv). For each such edge, 
Greedy ensures the existence of a (1 + e)-spanner path in Q. It follows that, at 
the end of iteration k, Q contains a (1 + e) fc -spanncr wi>-path. □ 

Theorem 7. [Leapfrog Property] Let L be an arbitrary grid cell and let F C 
El be the set of edges with both endpoints in L that get added to H in step 
3 of the algorithm. Then F satisfies the (1 + e,t)-leapfrog property, for t = 
(l + £) 4 (f + l)C del . 

Proof. Consider an arbitrary subset S — U2V2, ■ ■ ■ , u m v m } C F. To prove 

inequality (6) for S, it suffices to consider the case when mvi is a longest edge 
in S. Define S' — {v m ui}U {v^ui + \ \ 1 < I < s}. Since Ui and Vi lie in L for each 
i, all edges from 5" lie entirely in L. Let ab e S' be arbitrary. If \ab\ > |wii>i|, 
then inequality (6) trivially holds, so assume that \ab\ < \u\Vi\. Next we show 
that Q contains an a6-path of length no greater than t\ab\ at the time {ui,vi} 
gets queried. We distinguish two cases: 

(i) ab G YDel. In this case ab gets queried in step 3 prior to uii>i, meaning that 
Q contains a path Pq(a,b) of length |Pg(a, b)\ < (1 + £) 4 |afe|, at the time 
U\Vi gets queried (by Lem. 4). 

(ii) ab £ YDel. By Lem. 2, YDel contains a path PYDei(a, b) of length 

\P™ e i(a,b)\ < — ^\ab\ (7) 
(1 + e) 4 

that contains only edges shorter than ab. For each edge pq G -PYDei(o, b), Q 
contains a path Pq(p, q) of length \Pq(p, q) < (1 + e) 4 |p<7|, at the time u\Vi 
gets queried (by Lem. 4). Thus we have that 

\P Q (a,b)\= E \P Q (P,q)\<(l+e) A J2 \Pl\<t\ab\ (8) 

pqePwei(a,b) P9G-PyDei(a,6) 

This latter inequality follows from (7). 

For 1 < k < s, let Pe be a shortest Vfi^+i-path in Q, and let P m be a shortest 
VmWi-path in Q. By the arguments above, such paths exists in Q at the time U1V1 
gets queried, and their stretch factor does not exceed t. Then P = P\® u-ivi © 
Pi®uzv-$® ■ ■ .®P m is a path from m to vi in Q, and lo(P) is no greater than the 



right hand side of the leapfrog inequality (6). Furthermore, uj(P) > (l+£)|iii«i|, 
otherwise the edge u\V\ would not have been added to H (and Q) in step 3 of 
the algorithm. This concludes the proof. □ 

Theorem 8. The output H generated by PLDS has total weight 0(lo(M ST)) . 

Proof. The proof is very similar to the proof of Thm. 3 and uses the results of 
Lem. 3 and Thm. 7. □ 




Fig. 4. Valid ranges for \&p H (u, v)\ < (l + e)|twj| queries (step 3 of the PLDS algorithm), 
illustrated for e = 1/2: query range for edge uv (left), for edge ab (middle), and for the 
entire grid cell L (right). 



Lemma 5. For any e < 2, the shortest path query |spg(u, u)| < (1 + e)|uu| in 
step 3 of the PLDS algorithm involves only those grid cells incident to the cell L 
containing uv. 

Proof. For a fixed edge uv, the locus of all points z with the property that 
\uz\ + \zv\ < (1 +e)|wu| is a closed ellipse A with focal points u and v. Clearly, a 
point exterior to A cannot belong to a (1 + e)-spanner path p(u, v) from u to v, 
so it suffices to limit the search for p(u, v) to the interior of A. Fig. 4 (left and 
middle) shows the search domains for edges corresponding to one diagonal (uv) 
and one side (ab) of a grid cell. For any grid cell L, the union of L and the search 
ranges for the two diagonals and four sides of L covers the search domain for any 
edge that lies entirely in L (see Fig. 4 right). It can be easily verified that, for 
e < 2, the search domain for L fits in the union of L and its eight surrounding 
grid cells. □ 

Theorem 9. The PLDS algorithm can be implemented in 0(1) rounds of com- 
munication. 

Proof. Computing LDel in step 1 of the algorithm takes at most 4 communication 
rounds [21]. As shown in the proof of Thm. 3, computing the clique cover in 
step 1 takes at most 8 rounds of communication. Step 2 of the algorithm is 
restricted to cliques. A node u belongs to at most 4 cliques. For each such 
clique, u executes step 2 locally, on the neighborhood collected in step 1. In a 
few rounds of communication, each node u is also able to collect the information 
on the grid cells incident to the ones containing u. By Lem. 5, this information 
suffices to execute step 4 of the algorithm locally. □ 



3 Conclusions 



We present the first localized algorithm that produces, for any given QUDG 
G and any e > 0, a (f + e)-spanner for G of maximum degree 0(1) and total 
weight 0(u>(MST)), in O(l) rounds of communication. We also present the first 
localized algorithm that produces, for any given UDG G, a planar 0(l)-spanner 
for G of maximum degree O(l) and total weight 0{u{MST)), in 0(1) rounds of 
communication. Both algorithms require the use of a Global Positioning System 
(GPS), since each node uses its own coordinates and the coordinates of its neigh- 
bors to take local decisions. Our work leaves open the question of eliminating the 
GPS requirement without compromising the quality of the resulting spanners. 
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Abstract. We present a new efficient localized algorithm to construct, 
for any given quasi-unit disk graph G = (V, E) and any e > 0, a (1 + e)- 
spanner for G of maximum degree O(l) and total weight 0(ui(M ST)), 
where ui(MST) denotes the weight of a minimum spanning tree for V . We 
further show that similar localized techniques can be used to construct, 
for a given unit disk graph G = (V, E), a planar Cde!(l+£)(l+§)-spanner 
for G of maximum degree 0(1) and total weight 0(uj{M ST)). Here Cdei 
denotes the stretch factor of the unit Delaunay triangulation for V . Both 
constructions can be completed in O(l) communication rounds, and re- 
quire each node to know its own coordinates. 



1 Introduction 

For any fixed a, < a < 1, a graph G = (V,E) is an a- quasi unit disk graph 
(cv-QUDG) if there is an embedding of V in the Euclidean plane such that, for 
every vertex pair u, v S V, uv G E if \uv\ < a, and \uv\ ^ E if \uv\ > 1. The 
existence of edges with length in the range (a, 1] is specified by an adversary. If 
a = 1, G is called a unit disk graph (UDG). a-QUDGs have been proposed as 
models for ad-hoc wireless networks composed of homogeneous wireless nodes 
that communicate over a wireless medium without the aid of a fixed infrastruc- 
ture. Experimental studies show that the transmission range of a wireless node 
is not perfectly circular and exhibits a transitional region with highly unreliable 
links [34] (see for example Fig. la, in which the shaded region represents the 
actual transmission range). In addition, environmental conditions and physical 
obstructions adversely affect signal propagation and ultimately the transmission 
range of a wireless node. The parameter a in the a-QUDG model attempts to 
take into account such imperfections. 

Wireless nodes are often powered by batteries and have limited memory 
resources. These characteristics make it critical to compute and maintain, at 
each node, only a subset of neighbors that the node communicates with. This 
problem, referred to as topology control, seeks to adjust the transmission power 
at each node so as to maintain connectivity, reduce collisions and interference, 
and extend the battery lifetime and consequently the network lifetime. 
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Different topologies optimize different performance metrics. In this paper we 
focus on properties such as planarity, low weight, low degree, and the spanner 
property. Another important property is low interference [5, 15, 30], which we do 
not address in this paper. A graph is planar if no two edges cross each other (i.e, 
no two edges share a point other than an endpoint). Planarity is important to 
various memoryless routing algorithms [16, 4]. A graph is called low weight if its 
total edge length, defined as the sum of the lengths of all its edges, is within a 
constant factor of the total edge length of the Minimum Spanning Tree (MST). 
It was shown that the total energy consumed by sender nodes broadcasting along 
the edges of a MST is within a constant factor of the optimum [31]. Low degree 
(bounded above by a constant) at each node is also important for balancing 
out the communication overhead among the wireless nodes. If too many edges 
are eliminated from the original graph however, paths between pairs of nodes 
may become unacceptably long and offset the gain of a low degree. This renders 
necessary a stronger requirement, demanding that the reduced topology be a 
spanner. Intuitively, a structure is a spanner if it maintains short paths between 
pairs of nodes in support of fast message delivery and efficient routing. We define 
this formally below. 

Let G = (V, E) be a connected graph representing a wireless network. For 
any pair of nodes u, v € V, let sp G (u,w) denote a shortest path in G from u to 
v, and let |sp G (w, v)\ denote the length of this path. Let H C G be a connected 
subgraph of G. For fixed t > 1, H is called a t- spanner for G if, for all pairs of 
vertices u, v G V, \sp H (u,v)\ < t ■ |sp G (w,v)|. The value t is called the stretch 
factor of H . If t is constant, then H is called a length spanner, or simply a 
spanner. A triangulation of V is a Delaunay triangulation, denoted by Del(V), 
if the circumcirclc of each of its triangles is empty of nodes in V. 

Due to the limited resources and high mobility of the wireless nodes, it is 
important to efficiently construct and maintain a spanner in a localized manner. 
A localized algorithm is a distributed algorithm in which each node u selects all 
its incident edges based on the information from nodes within a constant number 
of hops from u. Our communication model is the standard synchronous message 
passing model, which ignores channel access and collision issues. In this commu- 
nication model, time is divided into rounds. In a round, a node is able to receive 
all messages sent in the previous round, execute local computations, and send 
messages to neighbors. We measure the communication cost of our algorithms in 
terms of rounds of communication. The length of messages exchanged between 
nodes is logarithmic in the number of nodes. 

Our Results. In this paper we present the first localized method to construct, 
for any QUDG G = (V, E) and any e > 0, a (1 + e)-spanncr for G of maximum 
degree 0(1) and total weight 0(uj{MST)), where uj{MST) denotes the weight 
of a minimum spanning tree for V. We further extend our method to construct, 
for any UDG G — (V,E), a planar spanner for G of maximum degree 0(1) and 
total weight 0(lu(MST)). The stretch factor of the spanner is bounded above by 
Cdez(L+ £ )(l. + § ), where Cdei is the stretch factor of the unit Delaunay triangula- 
tion for V (Cdei < 2.42 [20]). This second result resolves an open question posed 



by Li et al. in [22]. Both constructions can be completed in 0(1) communication 
rounds, and require each node to know its own coordinates. 



1.1 Related Work 

Several excellent surveys on spanners exist [27,26,14,25]. In this section we 
restrict our attention to localized methods for constructing spanners for a given 
graph G = (V,E). We proceed with a discussion on non-planar structures for 
UDGs first. Existing results are summarized in the first four rows of Table 1. 

The Yao graph [33] with an integer parameter k > 6, denoted YGk, is defined 
as follows. At each node u £ V, any k equal-separated rays originated at u define 
k cones. In each cone, pick a shortest edge uv, if there is any, and add to YGk the 
directed edge uv. Ties are broken arbitrarily or by smallest ID. The Yao graph is a 
spanner with stretch factor 1 _ 2s h 17r /fc > however its degree can be as high as n— 1. 
To overcome this shortcoming, Li et al. [18] proposed another structure called 
YaoYao graph YYk, which is constructed by applying a reverse Yao structure 
on YGk- at each node u in YGk, discard all directed edges vu from each cone 
centered at u, except for a shortest one (again, ties can be broken arbitrarily 
or by smallest ID). YYk has maximum node degree 2k, a constant. However, 
the tradeoff is unclear in that the question of whether YYk is a spanner or not 
remains open. Both YGk and YYk have total weight 0(n) ■ ui(MST) [6]. Li 
et al. [32] further proposed another sparse structure, called YaoSink YSk, that 
satisfies both the spanner and the bounded degree properties. The sink technique 
replaces each directed star in the Yao graph consisting of all links directed into a 
node u, by a tree T[u) with sink u of bounded degree. However, neither of these 
structures has low weight. 
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Table 1. Results on localized methods for UDGs. 



We now turn to discuss planar structures for UDGs. The relative neighbor- 
hood graph (RNG) [29] and the Gabriel graph (GG) [12] can both be constructed 
locally, however neither is a spanner [2]. On the other hand, the Delaunay tri- 
angulation Del(V) is a planar t-spanner of the complete Euclidean graph with 
vertex set V. This result was first proved by Dobkin, Friedman and Supowit [11], 
for t = 1+ 2 v ^ 7r w 5.08, and was further improved to t = ^p7r ~ 2.42 by Keil 
and Gutwin [17]. Das and Joseph [7] generalize these results by identifying two 
properties of planar graphs, the good polygon and diamond properties, which 
imply that the stretch factor is bounded above by a constant. 

For a given point set V, the unit Delaunay triangulation of V , denoted 
UDel(V), is the graph obtained by removing all Delaunay edges from Del(V) 
that are longer than one unit. It was shown that UDel(V) is a t-spanner of the 
unit-disk graph UDG0O, with t = ^tt w 2.42 [20]. 

Gao et al. [13] present a localized algorithm to build a planar spanner called 
restricted Delaunay graph (RDG), which is a supergraph of UDel(V). Li et al. [20] 
introduce the notion of a k-localized Delaunay triangle: Aabc is called k-localized 
Delaunay if the interior of its circumcircle does not contain any node in V 
that is a fc-neighbor of a, b or c, and all edges of Aabc are no longer than one 
unit. The authors describe a localized method to construct, for fixed k > 1, 
the k- localized Delaunay graph LDel k (V), which contains all Gabriel edges and 
edges of all k- localized Delaunay triangles. They show that (i) LDel k (V) is a 
supergraph of UDel(V) (and therefore a ^p7r-spanncr), (ii) LDel fe (V) is pla- 
nar, for any k > 2, and (iii) LDel 1 (F) may not be planar, but a planar subgraph 
PLDel(F) C LDel 1 (V) that retains the spanner property can be locally extracted 
from LDel 1 (T^). Their planar spanner constructions take 4 rounds of communi- 
cation and a total of 0(n) messages (O(nlogn) bits). Araujo and Rodrigues [1] 
improve upon the communication time for PLDel and devise a method to com- 
pute PLDel(V) in one single communication step. Both PLDel(T^) and LDel fc (y), 
for k > 1, may have arbitrarily large degree and weight. 

To bound the degree, several methods apply the ordered Yao structure on 
top of an unbounded-degree planar structure. This idea was first introduced by 
Bosc et al. in [3], and later refined by Li and Wang in [32, 23]. Since the ordered 
Yao structure is relevant to our work in this paper as well, we pause to discuss 
the Ordered Yao method for constructing this structure. The Ordered Yao 
method is outlined in Table 2. The main idea is to define an ordering ir of the 
nodes such that each node u has a limited number of neighbors (at most 5) who 
are predecessors in 7r; these predecessors are used to define a small number of 
open cones centered at u, each of which will contain at most one neighbor of u 
in the final structure. To maintain the spanner property of the original graph, 
a short path connecting all neighbors of u in each cone is used to replace the 
edges incident to u that get discarded from the original graph. 

Thm. 1 summarizes the important properties of the structure computed by the 
Ordered Yao method. 



Theorem 1. If G is a planar graph, then the output G' obtained by executing 
Ordered Yao(G) is a planar (1 + ■§)- spanner for G of maximum degree 25 [32]. 



Algorithm Ordered Yao(G = (V,E)) [32] 



{1. Find an order n for V:} 

Initialize i = 1 and G; = G. 
Repeat for i = 1, 2, . . . , \V\ 

Remove from G; the node u of smallest degree 

(break ties by smallest ID.) 

Call the remaining graph Gj+i. 

Set ty u = n — i + 1. 

{2. Construct a bounded-degree structure for G:} 

Mark all nodes in V unprocessed. Initialize E' <— and G' = (V, -E'). 
Repeat \V\ times 

Let u be the unprocessed node with the smallest order ir u . 
Let vi, V2, ■ ■ ■ , Vh be the be the processed neighbors of u in G (h < 5). 
Shoot rays from u through each Vi, to define /i sectors centered at u. 
Divide each sector into fewest open cones of degree at most it/3. 
For each such open cone C u (refer to Fig. above) 

Let si, S2, ■ ■ ■ , s m be the geometrically ordered neighbors of u in C u . 

Add to E' the shortest us; edge. 

Add to E' all edges SjSj+i, for j = 1, 2, . . . , m — 1. 
Mark node u processed. 

Output G = (V,E r ). 



Table 2. The Ordered Yao method. 




Song et al. [28] apply the ordered Yao structure on top of the Gabriel graph 
GG(V) to produce a planar bounded-degree structure OrdYaoGG. Their result 
improves upon the earlier localized structure YaoGG [18], which may not have 
bounded degree. Both YaoGG and OrdYaoGG are power spanners, however neither 
is a length spanner. 

The first efficient localized method to construct a bounded-degree planar 
spanner was proposed by Li and Wang in [32,23]. Their method applies the 
ordered Yao structure on top of LDel(V) to bound the node degree. The resulted 
structure, called BPS(V) (Bounded-Degree Planar Spanner), has degree bounded 
above by 19 + [^] , where < a < | is an adjustable parameter. The total 
communication complexity for constructing BPS(F) is 0(n) messages, however 
it may take as many as 0{n) rounds of communication for a node to find its rank 



in the ordering of V (a trivial example would be n nodes lined up in increasing 
order by their ID). The BPS structure does not have low weight [19]. 

The first localized low-weight planar structure was proposed in [19]. This 
structure, called RNG' , is based on a modified relative neighborhood graph, and 
satisfies the planarity, bounded-degree and bounded-weight properties. A similar 
result has been obtained by Li, Wang and Song [22], who propose a family of 
structures, called Localized Minimum Spanning Trees LMSTfc, for k > 1. The 
authors show that LMSTfc is planar, has maximum degree 6 and total weight 
within a constant factor of u>(MST), for k > 2. Their result extends an earlier 
result by Li, Hou and Sha [24], who propose a localized MST-based method to 
compute a local minimum spanning tree structure. However, neither of these 
low- weight structures satisfies the spanner property. Constructing low- weight, 
low-degree planar spanners in few rounds of communication is one of the open 
problems we resolve in this paper. 



2 Our Work 



We start with a few definitions and notation to be used through the rest of the 
paper. For any nodes u and v, let uv denote the edge with endpoints u and v; 
uv is the edge directed from u to v; and \uv\ denotes the Euclidean distance 
between u and v. Let C u denote an arbitrary cone with apex u, and let C u (v) 
denote the cone with apex u containing v. For any edge set E and any cone C u , 
let E fl C u denote the subset of edges in E incident to u that lie in C u . 

We assume that each node u has a unique identifier ID(tt) and knows its 
coordinates (x u ,y u )- Define the identifier ID(mu) of a directed edge uv to be 

the triplet (\uv\, ID(tt), ID(u)). For any pair of directed edges uv and u'v', we 
say that ID(uv) < ID(wV) if and only if one of the following conditions holds: 
(1) \uv\ < \u'v'\, or (2) \uv\ = \u'v'\ and ID(u) < ID(u'), or (3) \uv\ = \u'v'\ 
and ID(m) = ID(u') and ID(u) < ID(u'). For an undirected edge uv, define 
IT)(uv) — min{lD(wu), ID(vu)}. Note that according to this definition, each edge 
has a unique identifier. 

Let H — (V, Eh) be an arbitrary subgraph of G = (V, E). A subset L u <zV 
is an r-cluster in H with center u if, for any v e L u , \sp H (u, v)\ < r. A set 
of disjoint r-clusters {L Ul ,L U2 , . . .} form an r-cluster cover for V in H if they 
satisfy two properties: (i) for i ^ j, \sp H (ui,Uj)\ > r (the r-packing property), 
and (ii) the union UiL Ui covers V (the r-covering property). 

For any node subset U C V, let G[U] denote the subgraph of G induced by 
U. A set of node subsets V\, V2, . ■ ■ C V is a clique cover for V if the subgraph 
of G[Vi] is a clique for each i, and U^ =1 Fi = V. 

The aspect ratio of an edge set E is the ratio of the length of a longest edge 
in E to the length of a shortest edge in E. The aspect ratio of a graph is defined 
as the aspect ratio of its edge set. 



2.1 The LOS Algorithm 



In this section we describe an algorithm called LOS (Localized Optimal Spanner) 
that takes as input an a-QUDG G = (V, E), for fixed < a < 1, and a value 
e > 0, and computes a (1 +e)-spanner for G of maximum degree O(l) and total 
weight 0(uj(MST)). The main idea of our algorithm is to compute a particular 
clique cover Vi, V2, ■ ■ ■ for V, construct a (1 + e)-spanner for each G[Vi], then 
connect these smaller spanners into a (1 + e)-spanner for G using selected Yao 
edges. In the following we discuss the details of our algorithm. 
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Fig. 1. (a) The a-QUDG model (b) Constructing a clique cover for V (c) Clique 
ordering. 



Let < < -j= and < S < /3/4 be small constants to be fixed later. 
To compute a clique cover for V , we start by covering the plane with a grid of 
overlapping square cells of size (3 x /3, such that the distance between centers of 
adjacent cells is — 28. Note that any two adjacent cells define a small band of 
width S where they overlap. The reason for enforcing this overlap is to ensure 
that edges not entirely contained within a single grid cell are longer than 5, i.e., 
they cannot be arbitrarily small. We identify each grid cell by the coordinates 
(i,j) of its upper left corner. Any two vertices that lie within the same grid cell 
are no more than a distance apart and therefore are connected by an edge in 
G. This implies that the collection of vertices in each non-empty grid cell can 
be used to define a clique element of the clique cover. We call this particular 
clique cover a ((3 , S)- clique cover. Let Vi,V2, ■ ■ ■ be the elements of the ((3,5)- 
clique cover for V. Note that, since 5 < f3/4, a node u can belong to at most 
four subsets Vi. 

Our LOS method consists of 4 steps. First we construct, for each G[Vi], a 
(1 + e)-spanner of degree O(l) and weight 0(u(MST(Vi)). Various methods for 
constructing Hi exist - for instance, the well-known sequential greedy method 
produces a spanner with the desired properties [8]. Second, we use the Yao 
method to generate (1 + e) -spanner paths between longer edges that span differ- 
ent grid cells. Third, we apply the reverse Yao step to reduce the number of Yao 
edges incident to each node. Finally, we apply a filtering method to eliminate 
all but a constant number of edges incident to a grid cell. This fourth step is 
necessary to ensure that the output spanner has bounded weight. These steps 



Algorithm LOS(G = (V,E),e) 

{1. Compute a (1 + e)-spanner cover:} 

Fix < (3 < -z= and < <5 < (3/4. 

Compute a (/3, (5)-clique cover Vi, V2, . . . for V. 

For each i, compute a (1 + e)-spanner Hi for G[Vi] using the method from [8]. 
Initialize H = Uji/j. Let Eo = {uv £ E \ uv £ G[Vi] for any i}. 

{2. Apply Yao on E< } :} 

Let k be the smallest integer satisfying cos 2zl — sin 2zl > _ 
For each node w, divide the plane into k incident equal-size cones. 
Initialize Ey <— 0. 

For each cone C u such that i?o fl C u is non-empty 

Pick the edge uv G E PI C u of smallest ID and add uv to £V . 

{3. Apply reverse Yao on Ey:} 
Initialize Eyy <— _Ey. 

For each cone C u such that Ey n C u is non-empty 

Discard from _Ey all edges vu G -Ey fl C u , but the one of smallest ID. 

{4. Select connecting edges from Eyy:} 

Pick r such that r < g±lKl±£ii£g« fl-Binfl)-(g+i+ e ) > where fl = 2?r/ ^ 

Compute an r-cluster cover for V in ff . 

Let Ei C £yy contain all Yao edges connecting cluster centers. Add E\ to _ff. 
Output H = (V,E H ). 



Table 3. The LOS algorithm. 



are described in detail in Table 3. Note that the Yao and reverse Yao steps are 
restricted to edges in the set E whose aspect ratio is bounded above by 1/5. 
The next three theorems prove the main properties of the LOS algorithm. 

Theorem 2. The output H generated by LOSfG, e) is a (1 + e)- spanner for G. 

Proof. Let uv € E be arbitrary. If uv e G[Vi] for some i, then Hi C H contains 
a (1 + e)-spanner uw-path (since Hi is a (1 + e)-spanner for G[Vi]). Otherwise, 
uv e E . The proof that H contains a (1 + e)-spanner uw-path is by induction on 
the ID of edges in E . Let uv e E n be the edge with the smallest ID and assume 
without loss of generality that TD(uv) = TD{uv). Since ID(uv) is smallest, uv 
gets added to Ey in step 2, and it stays in Eyy in step 3. If uv <G H at the end 
of step 4, then sp H (u, v) — uv. Otherwise, let ab be the edge selected in step 4 
of the algorithm, such that u G L a and v € L& (see Fig. 2a). Since L a and Lb are 
both r-clusters, we have that \sp H (u, a)\ < r and \sj> H (v,b)\ < r. It follows that 
\ua\ < r and \vb\ < r. By the triangle inequality \ab\ < \uv\ + 2r and therefore 



sp ff (u, v) < \ab\ + 2r < \uv\ + Ar < (1 + e)|w|, for any r < Se/4 (satisfied by 
the r values restricted by the algorithm). This concludes the base case. 

To prove the inductive step, let uv € Eq be arbitrary, and assume that H 
contains (1 + e)-spanner paths between the endpoints of any edge whose ID is 
lower than ID(uv). 




(a) (b) 

Fig. 2. Thm. 2: (a) Base case, (b) sp H (it, ui) ©sp H (iti, a) ©afc©sp /J -(6, Vi) © sj> H (vi, v) 
is a (1 + e)-spanner wv-path. 



Let uvi G C u (v) be the Yao edge selected in step 2 of the algorithm; let 
uivi € C Vl (u) be the Yao Yao edge selected in step 3 of the algorithm; and let 
ab 6 H be the edge added to H in step 4 of the algorithm, such that ui € L a 
and vi € Lb (see Fig. 2b). Note that u and m may be disjoint or may coincide, 
and similarly for v and v±. In either case, the chains of inequalities ID(ui^i) < 
ID(wi) < IV (uv) and < < \uv\ hold. Let u[ be the projection of u\ 

on uv\. By the triangle inequality, 

\uu\\ < \uu[\ + \u[ui\ — \uvi\ — \uiVi\ + \uiUi\ < \uvi\ — \uivi\co$6+\uivi\siri6. 

(1) 

Similarly, if v[ is the projection of v\ on uv, we have 

\v\v\ < \vv[\ + \v[vx\ = \uv\ — \uv[\ + \v[vi\ < \uv\ — \uvi\ cos + \uvi\ sin9. (2) 

Since \uu\\ < \uvi\ < \uv\ and \viv\ < \uv\, by the inductive hypothesis H 
contains (l + e)-spanner paths sp H (u, Ui) and sp H (vi,v). Let Pj = sp H (u, Ui)(B 
sp H (vi, v). The length of Pi is 

\Pi\ <(l + e)-(\uu 1 \ + \v 1 v\). 

Substituting inequalities (1) and (2) yields 

\Pi\ < (l + £)\uv\ + (l + £)\uv 1 \(l-cos6 + sm0)-(l+e)\uiv 1 \(cos6-sm6). (3) 

Next we show that the path P = P\® sj> H (ui,a) (Bab® sp H (fe, V\) is a (1 + e)- 
spanner path from u to v in H, thus proving the inductive step. Using the fact 
that \ab\ < 2r + \uiV\\, \sj> H (u±, a)\ < r and \sp H (b, vi)\ < r, we get 



\P\ < |Pi| + K«x|+4r. 



(4) 



Substituting further > 5 and \uv\\ < 1 in (3) and (4) yields 

|P| < (1 + e)\uv\ + (4r + (1 + s)(l - cos 6 + sin 6) - 5(1 + s) (cos 6 - sin 6) - 5). 

Note that the second term on the right side of the inequality above is non-positive 
for any r and satisfying the conditions of the algorithm: 



Before proving the other two properties of H (bounded degree and bounded 
weight), we introduce an intermediate lemma. For fixed c > 0, call an edge 
set F c-isolated if, for each node u incident to an edge e e F, the closed disk 
disk^c) centered at u of radius c contains no other endpoints of edges in F. 
This definition is a variant of the isolation property introduced in [10]. Das et al. 
show that, if an edge set F satisfies the isolation property, then u)(F) is within 
a constant factor of the minimum spanning tree connecting the endpoints of F. 
Here we prove a similar result. 

Lemma 3. Let F be a c-isolated set of edges no longer than 1. Then u)(F) = 
0(1) ■ oj(T), where T is the minimum spanning tree connecting the endpoints of 
edges in F. 

Proof. Let P be a Hamiltonian path obtained by a taking a preorder traversal 
of T. If each edge uv e P gets associated a weight value u(uv) = |sp T (w, v)\, 
then it is well-known that uj(P) < 2uo(T). So in order to prove that w(F) is 
within a constant factor of w(T), it suffices to show that uj(F) = 0(ui(P)). Since 
F is c-isolated, the distance between any two vertices in T is greater than c and 
therefore w(P) > (n — l)c. On the other hand, no edge in F is greater than 1 
and therefore uo(F) < n. It follows that uu(F) = 0(w(P)). □ 

Theorem 4. The output H generated by running LOS fG, t) has maximum de- 
gree 0(1) and total weight 0(1) • lu(MST). 

Proof. The fact that H has maximum degree O(l) follows immediately from 
three observations: (a) each spanner Hi constructed in step 1 of the algorithm 
has degree O(l) [8], (b) a node u belongs to at most four subgraphs Hi, and (c) 
a node u is incident to a constant number of Yao edges (at most 2k) [18]. 

We now prove that the total weight for H is within a constant factor of 
ui(MST), which is optimal. The main idea is to partition the edge set Eh into 
a constant number of subsets, each of which has low weight. Consider first the 
(1 + e)-spanners constructed in step 1 of the algorithm. Each (1 + £r)-spanner Hi 
corresponds to a grid cell (i,j). Let F denote the set of edges in UeHe. Define the 
edge set F s C F to contain all spanner edges corresponding to those grid cells 
(i, j) whose indices i and j satisfy the condition (i mod 3) x 3 + j mod 3 = s. 
Intuitively, if two edges e\, e 2 € F s lie in different grid cells, then those grid cells 




This completes the proof. 



□ 



are separated by at least two other grid cells (see Fig. lc). This further implies 
that the closest endpoints of e\ and e-i are distance a or more apart. Also notice 
that it takes only 9 subsets F\, F2, . . . , Fg to cover F. 

Next we show that lj(F s ) = 0(u(T s )) for each s = 1, 2, . . . , 9, where T s is a 
minimum spanning tree connecting the endpoints in F s . To see this, first observe 
that F s combines the edges of several low-weight (1 + e)-spanners H Sl , H S2 , . . . 
with the property that u>(H Sl ) = 0(u>(T Sl )), where T 8l is a minimum spanning 
tree connecting the nodes in H Sl . Thus, in order to prove that ui(F s ) = 0(u(T 8 )), 
it suffices to show J2i = 0(w(T s )). We will in fact prove that 

£>(T Si ) < w(T a ) 

i 

We prove this by showing that, if Prim's algorithm is employed in constructing T s 
and T Si , then T Si C T s , for each i. Since the trees T Si are all disjoint (separated by 
at least 2 grid cells), the claim follows. Recall that Prim's algorithm processes 
edges by increasing length and adds them to T s as long as they do not close 
a cycle. This means that all edges shorter than a are processed before edges 
longer than a. Let e G T Si be arbitrary. Then |e| < a, since T Si is restricted 
to one grid cell only of diameter a. If e ^ T s , then it must be that e closes a 
cycle C at the time it gets processed. Note however that C must lie entirely in 
the grid cell containing T Si , since C contains edges no longer than a, and all 
edges with endpoints in different cells are longer than a. Furthermore, C must 
contain an edge e' ^ T Si such that |e'| < |e|. The case |e'| = \e\ cannot happen 
if Prim breaks ties in the same manner in both T s and T Si , so it must be that 
|e'| < \e\. But then we could replace e in T Si by e', resulting in a smaller spanning 
tree, a contradiction. It follows that e € T s and therefore T Si C T s , for each i. 
This concludes the proof that lj(F s ) = 0(oj(T s )) 7 for each s. Since there are 
at most 9 such sets F s that cover F and since oj(T s ) < uj(MST), we get that 
Lo(F) = 0(lo(MST)). 

It remains to prove that u{E H \ F) = 0{lo{MST)). Let d < 2k be the 
maximum number of edges in Eh \ F incident to any node in H. Partition 
the edge set Eh \ F into no more than 2d < Ak subsets Ei 7 E 2} ■ ■ ., such that 
no two edges in Ei share a vertex, for each i. We now show that co(Ei) = 
0(uj(MST)), for each i. Since there are only a constant number of sets Ei (4fc 
at most), it follows that uj{E h \ F) = 0(lj(MST)). The key observation to 
proving that iv(Ei) = 0{uj(MST)) is that any two edges uv, ab € E t have their 
closest endpoints - say, u and a - separated by a distance of at least r/t. This 
is because t\ua\ > \sp H (u, a)\ > r; the first part of this inequality follows from 
the spanner property of H, and the second part follows from the fact that u 
and a are centers of different r-clusters (a property ensured by step 4 of the 
algorithm). This implies that Ei is r/t- isolated, and by Lem. 3 we have that 
u{Ei) = 0(uj{MST)). 

We have established that u){F) = 0{uo{MST)) &nduo{E H \F) = 0(lu(MST)). 
It follows that w(H) = w(E H ) = 0(uj(MST)) and this completes the proof. □ 



Theorem 5. The LOS algorithm can be implemented in 0(1) rounds of com- 
munication using messages that are O(logn) bits each. 

Proof. Let x u and y u denote the coordinates of a node u. At the beginning of the 
algorithm, each node u broadcasts the information (lD(u), x u , y u ) to its neighbors 
and collects similar information from its neighbors. Each node u determines the 
grid cell(s) (i, j) it belongs to from two conditions, ia/^/2 < x u < (i + l)a/y/2 
and ja/\/2 < y u < (j + l)a/\/2. Similarly, for each neighbor v of u, each node 
u determines the grid cell(s) that v belongs to. Thus step 1 of the algorithm can 
be implemented in one round of communication: using the information from its 
neighbors, each node u computes the clique corresponding to those cells 
that u belongs to (at most 4 of them), then u computes a (1 + £)-spanner for 
each such clique by performing local computations. Note that knowledge of node 
coordinates is critical to implementing step 1 efficiently. 

Step 2 (the Yao step) and step 3 (the reverse Yao step) of the algorithm are 
inherently local: each node u computes its incident Yao and Yao Yao edges based 
on the information gathered from its neighbors in step 1. 

It remains to show that step 4 can also be implemented in 0(1) rounds of 
communication. We will in fact show that eight rounds of communication suffice 
to compute an r-cluster cover for V in H . Define U s to be the set of vertices 
that lie in the grid cells such that (i mod 2) x 2 + j mod 2 = s. This is 
the same as saying that two vertices that lie in different cells are about one grid 
cell apart. Note that V — Uf =1 ?7 s . To compute an r-cluster cover for V, each 
node u executes the ClusterCover method described below. For simplicity we 
assume that r > 5, so that two cluster centers that lie in different grid cells are 
at least distance r apart. However, the ClusterCover method can be easily 
extended to handle the situation r < S as well. 



Computing a ClusterCover^, r) 

Repeat for s = 1, 2, 3, 4 

(A) Collect information on cluster centers from neighbors (if any). 
If u belongs to U s 

Let Ve C Us be the clique containing u (computed in step 1 of LOS). 

(B) Broadcast information on existing cluster centers in Ve to all nodes in V(. 

(C) For each existing cluster center w € Ve 

Add to C m all uncovered nodes v £ Ve such that sp H (w,v) < r. 
Mark all nodes in C w covered. 

(D) While Ve contains uncovered nodes 

Pick the uncovered node w € Ve of highest ID. 

Add to C w all uncovered nodes v £ Ve such that sTp H (v,w) < r. 

Mark all nodes in C w covered. 

(E) Broadcast the cluster centers computed in step (C) to all neighbors. 



No information on existing cluster centers is available in the first iteration of the 
ClusterCover method (i.e, for s = 1). Each node in Ui skips directly to step 



(D), which implements the standard greedy method for computed an r-clique 
cover for a given node set (Vg in our case). In the second iteration, some of the 
clusters computed during the first iteration might be able to grow to incorporate 
new vertices from L^- This is particularly true for cluster centers that lie in the 
overlap area of two neighboring cells. Information on such cluster centers is 
distributed to all relevant nodes in step (E) in the first iteration, then collected 
in step (A) and forwarded to all nodes in Vi in step (B) in the second iteration. 
This guarantees that all nodes in have a consistent view of existing cluster 
centers in Ve at the beginning of step (C). Existing clusters grow in step (C), if 
possible, and new clusters get created in step (D), if necessary. This procedure 
shows that it takes no more than 8 rounds of communication to implement step 
4 of the LOS algorithm. One final note is that information on a constant number 
of cluster centers is communicated among neighbors in steps (A), (B) and (D) 
of the ClusterCover method. This is because only a constant number of r- 
clusters can be packed into a grid cell. So each message is O(logn) bits long, 
necessarily so to include a constant number of node identifiers, each of which 
takes O(logn) bits. □ 



2.2 The PLOS Algorithm 

In this section we impose our spanner to be planar, at the expense of a bigger 
stretch factor. This tradeoff is unavoidable, since there are UDGs that contain 
no (1 + e)-spanner planar subgraphs, for arbitrarily small e (a simple example 
would be a square of unit diameter). 

Our PLOS algorithm consists of 4 steps. In a first step we construct the unit 
Delaunay triangulation UDel(V) using the method described in [21]. Remaining 
steps use the grid-based idea from Sec. 2.1 to refine the Delaunay structure. 
Let Vi,V2,... be a (/3, <S)-clique cover for V, as defined in Sec. 2.1. In step 2 
of the algorithm we apply the Ordered Yao method on edge subsets of UDel 
incident to each clique Vi . The reason for restricting this method to each clique, 
as opposed to the entire spanner Wel(V) as in [32], is to reduce the total of 0(n) 
rounds of communication to O(l). The individual degree of each node increases 
as a result of this alteration, however it remains bounded above by a constant. 
Steps 3 and 4 aim to reduce the total weight of the spanner. Step 3 uses a Greedy 
method to filer out edges with both endpoints in one same clique Vi . Step 4 uses 
clustering to filter out edges spanning multiple cliques. These steps are described 
in detail in Table 4. The reason for breaking up step 3 of the algorithm into 4 
different rounds (for k = 1, ... ,4) will become clear later, in our discussion of 
communication complexity (Thm. 15). We now turn to proving some important 
properties of the output spanner. We start with a preliminary lemma. 

Lemma 6. The graph YDel constructed in step 2 of the PLOS algorithm is a 
planar t\- spanner for G, for any t\ > G'deii^ + 1)- Furthermore, for each edge 
ah e G, YDel contains a t\-spanner ab-path with all edges shorter than ab [32]. 



Algorithm PLDS(G = (V, E),e) 

{1. Start with the localized Delaunay structure for G:} 
Compute LDel = (V,E LBel ) for G using the method from [21]. 
Fix < (3 < -±= and < <5 < f . Compute a (/?, <5)-clique cover Vi, Vi, . . . for V. 

{2. Bound the degree:} 

For each clique Vi do the following: 

2.1 Let Ei C -BuDei contain all unit Delaunay edges incident to nodes in Vi. 

2.2 Execute YDeli <- Ordered YAo(Gi = (V,Ei)) (see Table 2). 
Set YDel = (V, E Wel ) = \J f YDel,. 

{3. Bound the weight of edges confined to single grid cells:} 

Initialize Eh = and H = (V, Eh). 
Repeat for k = 1, 2, 3,4 

{Use Greedy on non-adjacent grid cells:} 

For each grid cell L = L(i,j) such that (i mod 2) x 2 + j mod 2 = k 

3.1 Let = £rDei H L contain all edges in YDel that lie in L. 

Let Eq — E mel \ El and Q — (V, Eq) define the query graph for El- 

3.2 Sort El in increasing order by edge ID. 

For each edge e = uv £ El, resolve a shortest path query: 
If spg(u,w) > (1 +e)|wu| then add uv to H and Q. 
Otherwise, eliminate uv from YDel. 

{4. Bound the weight of edges spanning multiple grid cells:} 

Pick r such that r < ^ and compute an r-cluster cover for YDel. 
Add to H those edges in YDel connecting cluster centers. 

Output H = (V,E H ). 



Table 4. The PLOS algorithm. 



Proof. LDel is a planar C^-spanncr for G [21]. By Thm. 1, YDel^ is a planar 
(| + l)-spanner for Gi, for each i. These together with the fact that LDel = \J t d 
show that YDel is a ii-spanner for G. 




(a) c w (b) 



Fig. 3. YDel is planar: edges ab and uv cannot cross. 

The fact that YDel is planar follows an observation in [32] stating that, if 
a non-Delaunay edge e e YDel crosses a Delaunay edge e', then e' must be 
longer than one unit and does not belong to YDel. More precisely the following 
properties hold: 

(a) A non-Delaunay edge ab € YDel cannot cross a Delaunay edge uv € YDel. 
Recall that each non-Delaunay edge ab € YDel closes an empty triangle Aabc 
whose other two edges ac and be are Delaunay edges. Thus, if ab crosses uv, 
then at least one of ac and be must cross uv, contradicting the planarity of 
LDel (see Fig 3a). 

(a) No two non-Delaunay edges ab,uv G YDel cross. The arguments here are 
similar to the ones above: if ab and uv intersect, then at least two of the 
incident Delaunay edges intersect, contradicting the planarity of LDel (see 
Fig. 3b). 

The second part of the lemma follows from [32]. □ 

Theorem 7. The output H generated by PLOS (G,e) is a planar t- spanner for 
G, for any constant t > Cd e ;(l + e)(l + §). 

Proof. Since H C YDel, by Lem. 6 we have that H is planar. We now show that 
H is a i-spanner for G. The proof is by induction on the length of edges in H . 
The base case corresponds to the edge uv e G of smallest ID. Clearly uv e LDel, 
since uv is a Gabriel edge. Also uv G YDel, since it has the smallest ID among all 
edges and therefore it belongs to the Yao structure for LDel. We now distinguish 
two cases: 

(a) There is a grid cell containing both u and v. In this case uv G H , since uv is 
the first edge queried by Greedy in step 3 and therefore it gets added to H. 

(b) There is no grid cell containing both u and v. Let ab be the edge selected 
in step 4 of the algorithm, such that m£L„ and v € Lb (see Fig. 2a). Then 
arguments similar to the ones used for the base case of Thm. 2 show that 
sp H (u, a) (Bab® sp H (b,v) is a (1 + e)-spanner uv-p&th, for any r < £(5/4. 



This concludes the base case. To prove the inductive step, pick an arbitrary edge 
uv G G, and assume that H contains i-spanner paths between the endpoints of 
each edge in G of smaller ID. By Lem. 6, YDel contains a y^-spanner path 
u = w , Mi, ...,u s = v: 

s t 

< rqr^lH ( 5 ) 

i=0 

For each edge UiU i+ i e YDel, one of the following cases applies: 

(a) There is a grid cell containing both m and Ui+\. In this case, the Greedy step 
(step 3 of the algorithm) guarantees that |sp ff (m, < (1 + e)\uiUi + i\. 

(b) There is no grid cell containing both Ui and u i+ \. Arguments similar to the 
ones for the base case show that |sp ff (m, Ui + i)\ < (1 + e)\uiU i+ i\. 

In either case, H contains a (1 + e)-spanner UiU i+ i-path. This together with (5) 
shows that 

s s 

|sp H (u, v)\ = ^ \ s Ph( u ^ u *+i)\ < (1 + e) X] ^ 

i=0 i=0 

This completes the proof. □ 

Theorem 8. T/ie output H generated by PLDS /ias maximum degree 0(1). 

Proof. Since H C YDel, it suffices to show that the graph YDel constructed 
in step 2 of the PLOS algorithm has degree bounded above by a constant. By 
Thm. 1, the maximum degree of YDel^ is 25, for each %. Also note that unit disk 
centered at a node u intersects 0(^2) grid cells, meaning that u is a neighbor of 
nodes in O(-^j-) grid cells and therefore belongs to a constant number of graphs 
YDel,;. This implies that the maximum degree of u is 25 • O(^j-), which is a 
constant. □ 

Definition 9. [Leapfrog Property] For any t > t' > 1, a set F of edges has 
the (t', t)-leapfrog property if, for every subset S = {uiVi,u 2 V2, ■ ■ ■ , u m v m } of F, 

t' ■ \uivi\ < X \uiVi\ +t ■ ( + \v m ui\y (6) 

j=2 i=l 

Das and Narasimhan [9] show the following connection between the leapfrog 
property and the weight of the spanner. 

Lemma 10. Let t > t' > 1. If the line segments F in d-dimensional space 
satisfy the (t' ,t) -leapfrog property, then lu(F) = 0{u>(MST)), where MST is a 
minimum spanning tree connecting the endpoints of line segments in F. 

Lemma 11. At the end of each iteration k in step 3 of the PLOS algorithm, for 
k = 1, . . . ,4, Q contains (1 + e) k -spanner paths between the endpoints of any 
YDel edge processed in iterations 1 through k. 



Proof. The proof is by induction on k. The base case corresponds to k = 1. In 
this case, Greedy ensures that Q contains a (l + e)-spanner uw-path for each edge 
uv processed in this iteration. This is because uv € YDel either gets added to Q 
in step 3.1 (and never removed thereafter), or gets queried in step 3.2. To prove 
the inductive step, consider a particular iteration k > 1, and assume that the 
lemma holds for iterations I = 1 . . . k — 1. Again Greedy ensures that Q contains 
a (l+e)-spanner iru-path for each edge uv processed in iteration k. Consider now 
an arbitrary edge uv processed in iteration I < k. By the inductive hypothesis, 
at the end of round k— 1, Q contains a (1 + e) fe_1 -spanner path p(u, v). However, 
it is possible that p(u, v) contains edges processed in round k (since Greedy does 
not restrict p{u, v) to lie entirely in the cell containing uv). For each such edge, 
Greedy ensures the existence of a (1 + e)-spanner path in Q. It follows that, at 
the end of iteration k, Q contains a (1 + e) fc -spanncr uv-path. □ 

Theorem 12. [Leapfrog Property] Let L be an arbitrary grid cell and let 
F C El be the set of edges with both endpoints in L that get added to H in 
step 3 of the algorithm. Then F satisfies the (1 + e,t)-leapfrog property, for 
t=(l + e) 4 (§ +l)C del . 

Proof. Consider an arbitrary subset S — {uiV\,U2V2, ■ ■ ■ , u m v m } C F. To prove 
inequality (6) for 5, it suffices to consider the case when u\V\ is a longest edge 
in S. Define S' — {v m ui}U {v^ui + \ \ 1 < I < s}. Since Ui and Vi lie in L for each 
i, all edges from 5" lie entirely in L. Let ab e S' be arbitrary. If \ab\ > |uii>i|, 
then inequality (6) trivially holds, so assume that \ab\ < |tii«i|. Next we show 
that Q contains an a6-path of length no greater than t\ab\ at the time {u\,vi} 
gets queried. We distinguish two cases: 

(i) ab G YDel. In this case ab gets queried in step 3 prior to UiV\, meaning that 
Q contains a path Pq(a,b) of length \Pq(cl, b)\ < (1 + £) 4 |afe|, at the time 
U\Vi gets queried (by Lem. 11). 

(ii) ab £ YDel. By Lem. 6, YDel contains a path PYDei(a, b) of length 

\P™ e i(a,b)\ < — ^\ab\ (7) 
(1 + e) 4 

that contains only edges shorter than ab. For each edge pq e PyDei(a, b), Q 
contains a path Pq(p, q) of length \Pq(p, q) < (1 + e) 4 |p<7|, at the time U\V\ 
gets queried (by Lem. 11). Thus we have that 

\P Q (a,b)\= E \P Q (P,q)\<(l+e) 4 J2 \Pl\<t\ab\ (8) 

pqePwei(a,b) P9G-PyDei(a,6) 

This latter inequality follows from (7). 

For 1 < k < s, let Pi be a shortest Vfi^+i-path in Q, and let P m be a shortest 
VmWi-path in Q. By the arguments above, such paths exists in Q at the time UiVi 
gets queried, and their stretch factor does not exceed t. Then P = P\® UiVi © 
-P2©«3^3© ■ ■ -©Pm is a path from u\ to v\ in Q, and lu(P) is no greater than the 



right hand side of the leapfrog inequality (6). Furthermore, uj(P) > (l+£)|iii«i|, 
otherwise the edge u\V\ would not have been added to H (and Q) in step 3 of 
the algorithm. This concludes the proof. □ 

Theorem 13. The output H generated by PLDS has total weight 0(u>(MST)). 

Proof. The proof is very similar to the proof of Thm. 4 and uses the results of 
Lem. 10 and Thm. 12. □ 




Fig. 4. Valid ranges for \sj> H (u, v)\ < (l + e)|twj| queries (step 3 of the PLDS algorithm), 
illustrated for e = 1/2: query range for edge uv (left), for edge ab (middle), and for the 
entire grid cell L (right). 



Lemma 14. For any e < 2, the shortest path query \sPq(u,v)\ < (l+e)|m;| in 
step 3 of the PLDS algorithm involves only those grid cells incident to the cell L 
containing uv. 

Proof. For a fixed edge uv, the locus of all points z with the property that 
\uz\ + \zv\ < (1 +e)|wu| is a closed ellipse A with focal points u and v. Clearly, a 
point exterior to A cannot belong to a (1 + e)-spanner path p(u, v) from u to v, 
so it suffices to limit the search for p(u, v) to the interior of A. Fig. 4 (left and 
middle) shows the search domains for edges corresponding to one diagonal (uv) 
and one side (ab) of a grid cell. For any grid cell L, the union of L and the search 
ranges for the two diagonals and four sides of L covers the search domain for any 
edge that lies entirely in L (see Fig. 4 right). It can be easily verified that, for 
e < 2, the search domain for L fits in the union of L and its eight surrounding 
grid cells. □ 

Theorem 15. The PLOS algorithm can be implemented in 0(1) rounds of com- 
munication. 

Proof. Computing LDel in step 1 of the algorithm takes at most 4 communication 
rounds [21]. As shown in the proof of Thm. 4, computing the clique cover in 
step 1 takes at most 8 rounds of communication. Step 2 of the algorithm is 
restricted to cliques. A node u belongs to at most 4 cliques. For each such 
clique, u executes step 2 locally, on the neighborhood collected in step 1. In a 
few rounds of communication, each node u is also able to collect the information 
on the grid cells incident to the ones containing u. By Lem. 14, this information 
suffices to execute step 4 of the algorithm locally. □ 



3 Conclusions 



We present the first localized algorithm that produces, for any given QUDG 
G and any e > 0, a (f + e)-spanner for G of maximum degree 0(1) and total 
weight 0(u>(MST)), in O(l) rounds of communication. We also present the first 
localized algorithm that produces, for any given UDG G, a planar 0(l)-spanner 
for G of maximum degree O(l) and total weight 0{u{MST)), in 0(1) rounds of 
communication. Both algorithms require the use of a Global Positioning System 
(GPS), since each node uses its own coordinates and the coordinates of its neigh- 
bors to take local decisions. Our work leaves open the question of eliminating the 
GPS requirement without compromising the quality of the resulting spanners. 
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